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Letter to the Editor 

The distance to the Pleiades from orbital solution of the 
double-lined eclipsing binary HD 23642^'** 



o 
o 

(N 



oo 



> 

m 
o 

o 
6 



U. Munari^, S. Dallaporta^, A. Siviero^'^, C. Soubiran'^, M. Fiorucci^, and P. Girard^ 

^ INAF Osservatorio Astronomico di Padova, Sede di Asiago, 1-36012 Asiago (VI), Italy 

^ Via Filzi 9, 1-38034 Cembra (TN), Italy 

^ Astrophysical Observatory of the Dept. of Astronomy, Univ. of Padova, 1-36012 Asiago (VI), Italy 

* Observatoire Aquitain des Sciences de I'Univers, CNRS UMR 5804, BP 89, 33270 Floirac, France 

Received date ; accepted date 



Abstract. Combining precise B,V photometry and radial velocities, we have been able to derive a firm orbital 
solution and accurate physical parameters for the newly discovered eclipsing binary HD 23642 in the Pleiades open 
cluster. The resulting distance to the binary and therefore to the cluster is 132±2 pc. This closely confirms the 
distance modulus obtained by classical main sequence fitting methods (m — M=5.60 or 132 pc), moving cluster 
techniques and the astrometric orbit of Atlas. This is the first time the distance to a member of the Pleiades 
is derived by orbital solution of a double-lined eclipsing binary, and it is intended to contribute to the ongoing 
discussion about the discordant Hipparcos distance to the cluster. 
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1. Introduction 

The distance to the Pleiades has been derived by main- 
sequence fitting methods several times over the years, and 
consistently found to cluster around m — Af =5.60 ±0.04 
(e.g. Turner 1979, Meynet et al. 1993), corresponding to a 
distance of 132 (±2) pc. Values of to - M=5.52 (127 pc) 
by Mitchell & Johnson (1957) and m - M=5.75 (141 pc) 
by Eggen (1950) bracket the range of published distance 
moduli. Therefore, it came as a surprise the shorter dis- 
tance of 116 ±3.3 pc (to — A/=5.33 ±0.06) derived from 
the parallaxes of 54 Pleiades members when the results of 
the Hipparcos astrometric mission became available (van 
Leeuwen & Hansen Ruiz 1997). The ~0.3 mag difference 
has far reaching consequences in many areas of astro- 
physics, and it soon prompted extensive observational and 
theoretical work to account for it. 
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* The photometric and spectroscopic observations discussed 
in this paper are electronically available from the web page 
jhttp : / /uUsse .pd. astro .it /HD _23642/ 

** Based in part on spectra collected with the ELODIE spec- 
trograph at the 1.93m telescope of the Observatoire de Haute 
Provence (OHP), France 



An anomalous abundance of Helium was discussed 
as a possible explanation of the large difference between 
Hipparcos and ground-based photometric distances to the 
Pleiades (Mermilliod ct al. 1997a). The required helium 
over-abundance (y=0.37) is however too large to be a fea- 
sible explanation according to Pinsonneault et al. (1998). 
A sub-solar metallicity could reconcile the distances, but 
the metallicity of the Pleiades has been measured several 
times and with different methods, which generally con- 
verged to a mean value of [Fe/H]=0.00±0.03 (e.g. Stauffer 
et al. 2003). Castellani et al. (2002) were able to obtain 
a good fit of the Pleiades main sequence with their the- 
oretical isochrones (with updated physical inputs) and 
the Hipparcos distance, employing a sub-solar metallicity 
of [Fe/H]=-0.15. Working differentially with a sample of 
nearby G and K stars, Percival et al. (2003) presented a re- 
calibration of photometric metallicity indexes and argued 
that they support a sub-solar metallicity of the Pleiades 
in spite of the spectroscopic evidences (e.g. Boesgaard & 
Friel 1990) for a solar one, therefore removing the dif- 
ference between main-sequence fitting and Hipparcos dis- 
tances. However, Stauffer et al. (2003) remarked about 
the anomalous blue colors of K stars in the Pleiades, that 
they ascribed to fast rotation and spotted surfaces as a 
consequence of the young age of the cluster. 
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orbital phase 



Fig. 1. Our photo-electric V (top) and B (bottom) observations of HD 23642 are phase plotted in the upper panel 
following the cphemeris of Table 2. The lower panel displays our radial velocities from Table 1 (filled circles) and those 
from Pearce (1957, squares) and Abt (1958, crosses). The curves over plot the orbital solution given in Table 2. 



Narayanan & Gould (1999) derived a distance modu- 
lus to the Pleiades of to — M=5.58 (±0.18) using a variant 
of the moving cluster method, the gradient in the radial 
velocity of the cluster members in the direction of the 
proper motion of the cluster. Their techniques relies on the 
assumption that the velocity structure of the Pleiades is 
not significantly affected by rotation. Narayanan & Gould 
concluded that the errors in the Hipparcos parallaxes to- 
ward Pleiades are spatially correlated over angular scales 
of a few degrees, with an amplitude up to 2 mas. Almost 
simultaneously, however, van Leeuwen (1999) compared 
the results on 9 clusters and concluded instead that the 
Hipparcos parallaxes of the Pleiades were basically unaf- 
fected by systematics and that the distance modulus to 
the latter is m - Af=5.37 (±0.07). Makarov (2002) sug- 
gested that, while Hipparcos parallaxes are overall of ex- 
cellent quality, those for the Pleiades suffered from data 
analysis procedures that could have introduced system- 
atics when a rich and bright cluster was crossing one 
of the two Hipparcos fields of view while the other was 
essentially deprived of stars, as in the case of Pleiades. 
From intermediate Hipparcos data and under some as- 
sumptions, Makarov has recomputed the distance to the 
Pleiades as 129 ±3 pc. F. van Leeuwen (private communi- 
cation) is currently working on a further iteration on orig- 
inal Hipparcos reductions that should assess and quantify 
is any systematics affect the Pleiades data. 

The question about the distance to the Pleiades has 
not yet been firmly solved, and to tackle it new, robust 
and independent approaches (as much geometric as pos- 
sible) are required. A first one, advocated by Paczynski 



(2003), combines astrometric and spectroscopic observa- 
tions of Atlas (HD 23850), one of the brightest members of 
the Pleiades, which is an astrometric binary with a period 
of 291 day, a semi-major axis of 12.9 mas and 0.246 eccen- 
tricity. Pan ct al. (2004) did not have the radial velocities, 
but they anyway derived the distance by combining the 
astrometric orbit with a mass-luminosity relation. They 
derived a distance of 135±2 pc, which is in close agree- 
ment with the results of main-sequence fitting methods. 

Another method involves double-lined eclipsing bina- 
ries (SB2 EB), which are the distance indicators now pro- 
viding the most reliable distances to Magellanic Clouds 
and other galaxies in the Local Group. The recent dis- 
covery by Torres (2003) that HD 23642 in the Pleiades 
is an SB2 eclipsing binary (the only one so far known in 
the cluster) prompted us to use it to measure the distance 
to the Pleiades. In this Letter we present accurate new B 
and V photoelectric photometry and radial velocities of 
HD 23642, and use them to derive an orbital solution and 
infer the distance to the star, and therefore to the cluster. 

2. Observations 

We observed HD 23642 in B and V (standard Johnson 
filters) from a private observatory near Cembra (Trento), 
Italy. The instrument was a 28 cm Schmidt-Cassegrain 
telescope equipped with an Optec SSP5 photometer. It 
proved already to be a very accurate and reliable in- 
strument (cf. Siviero et al. 2004) perfectly suited to 
deal with the low amplitude of HD 23642 eclipses (~0.1 
mag). HD 23568 (HIP 17664, Ft=6.824±0.011, Bt = 
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Table 1. Heliocentric radial velocities of HD 23642. The 
S/N is evaluated around 5500 A. 



date UT HJD 

2004 (+2453000) 

Feb 3 21.8689 39.41273 

Feb 3 22.8489 39.45356 

Feb 4 0.0949 39.50548 

Feb 7 17.9842 43.25049 

Feb 10 18.0764 46.25318 



S/N RVi RVs 

(km s-i) (km s-i) 

152 -85.0±0.5 +134.5±0.5 

118 -88.0±0.5 +139.0±0.5 

84 -91.0±0.5 +143.0±1.0 

106 +105.5±1.5 -136.0±1.0 

121 +29.5±0.5 -24.5±1.0 



6.842 ±0.015, spectrum B9.5V) was chosen as comparison 
star and HD 23763 (HIP 17791, yT=6.963±0.011, Bt = 
7.109 ± 0.016, spectrum AlV) as a check star. Following 
the Bessell (2000) transformations between Tycho and 
Johnson photometric systems, we adopted F=6.830 and 
5=6.851 for the comparison star. The comparison star 
was measured against the check star at least once every 
observing run, and found constant with standard devia- 
tions of 0.005 mag in V and 0.006 mag in B, confirming 
the Hipparcos photometric results. In all, 492 measure- 
ments in and 432 in B were collected of HD 23642 
between Aug. 19, 2003 and Feb. 16, 2004. All the obser- 
vations were corrected for atmospheric extinction and in- 
strumental color equations (via calibration on Landolt's 
equatorial fields), and the instrumental differential mag- 
nitudes were transformed into the standard Johnson UBV 
system. The variable, comparison and check stars are sim- 
ilar in spectral type, lie very close on the sky (15 arcmin) 
and were always observed at zenith distances < 60°, which 
argues for a high consistency of our photometry. 

High resolution spectra of HD 23642 were obtained at 
five distinct epochs with the ELODIE echelle spectrograph 
at the 1.93m telescope of the Haute-Provence Observatory 
(Baranne et al. 1996). ELODIE covers the spectral range 
3900-6800 A in a single exposure as 67 orders at a mean 
resolving power of 42000. Optimal extraction and wave- 
length calibration were performed with the on-line au- 
tomatic reduction pipeline. Radial velocities for the two 
components are reported in Table 1. They were measured 
by both classical cross-correlation against suitable A-type 
templates and by the Least-Square Deconvolution method 
of Donati et al (1997). The errors given in Table 1 reflect 
the scatter of the measurements obtained with the two 
methods run with different parameters. 

3. The temperature of the primary and the 
reddening 

One parameter, the temperature of the primary (Ti), can- 
not be directly modeled in SB2 EB analysis and must be 
determined independently. In the case of HD 23642, the 
spectral classification reported in literature is too sparse, 
ranging from B9 to Al. To determine Ti we then turned to 
photometry and searched the literature and the General 
Catalogue of Photometric Data (GCPD, Mermilliod et 



Table 2. Orbital solution for HD 23642 (over plotted to 
observed data in Figure 1). Formal errors of the solution 
are given, f: derived from analysis of multi-band photom- 
etry and fixed. 
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al. 1997b) for observations of HD 23642. Data were 
foimd in the Walraven VBLUW, Stromgren ubvy, WBVR, 
Geneva, DDO, Vilnius, Johnson UBVRL Tycho, 2MASS 
and ASN 74 photometric systems (Crawford & Perry 1976, 
Kornilov et al. 1991, Morel & Magnenat 1978, Persinger 
& Castelaz 1990, van Lceuwen et al. 1986, Wesselius et 
al. 1982, GCPD). Using the relevant parameters for these 
systems as censed in the Asiago Database on Photometric 
Systems (Moro & Munari 2000, Fiorucci & Munari 2003), 
they were applied to the extensive synthetic spectral atlas 
of Munari et al. (2004). A total of 2310 synthetic spectra 
of the binary were retrieved from it considering all possi- 
ble combinations of 7 different temperatures for the pri- 
mary (9000< Ta < 10500 K), 10 for the secondary (6750< 
Ta <9000 K), 3 values for the gravity (log. 9^=4.0, 4.5, 
5.0) and 11 for the reddening (0.00< Eb-v <0.10). The 
[Fe/H]=0.00±0.02 metalhcity of the Pleiades reported by 
Stauffer et al. (2003) was adopted, as well as the stan- 
dard Rv = Av/EB-v=-i-^ reddening law (Fitzpatrick 
1999) that was found to be valid for the Pleiades by 
Guthrie (1987). A minimum distance method was then 
applied to observed and computed colors in each photo- 
metric system, which provided pretty consistent results 
among them. The weight averaged means for the temper- 
ature of the primary and the reddening were well con- 
strained to Ti = 9671 ± 46 K and £;B^y=0.012±0.004. 
Less constrained were the temperature of the fainter sec- 
ondary (T2 = 8023 ± 544 K) and the surface gravity of 
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the primary (log (?i=4.17±0.13), while the gravity of the 
secondary gave no measurable signal. 

The value of the reddening is perfectly confirmed by 
the equivalent width of the interstellar Nal D lines. We 
measured them on our high resolution spectra and found 
E.W.(NaI 5890 A)=0.035±0.001 A that corresponds to 
i;B_y=0.011 adopting the Munari & Zwitter (1997) cal- 
ibration. White et al. (2001) has derived a closely similar 
equivalent width (0.036±0.003 A) from a single, higher 
resolution spectrum. 

4. Orbital solution and distance to the Pleiades 

To augment the number and phase coverage of the ra- 
dial velocities, in the orbital solution we also considered 
the measurements by Pearce (1957) and Abt (1958). The 
adopted average relative weights are 10, 2, 1.5 for ours, 
Pearce and Abt data, respectively. 

The orbital solution of HD 23642 has been obtained 
with version WD98K93d of the Wilson-Devinney code 
(Wilson & Devinney 1971) as modified by Milone et al. 

(1992) to include Kurucz's model atmospheres to approx- 
imate the surface fluxes of the two stars. The computa- 
tions have been run within Mode-2 program option. Limb 
darkening coefficients have been taken from Van Hamme 

(1993) interpolated for the metallicity, temperature and 
gravity appropriate for the components of HD 23642. A 
linear law for limb darkening and single reflection (i.e. 
only the inverse square law illumination) were considered. 
The bolometric albedos were set to 1.0 and tests on the 
lightcurve confirmed the choice. Finally, a gravity bright- 
ening exponent f3— 1.0 has been adopted consistent with 
the radiative nature of the atmospheres in HD 23642. The 
effects on the orbital solution of the adopted limb dark- 
ening coefficients and law, bolometric albedos, single re- 
flection scenario and gravity brightening exponent are es- 
sentially negligible. Their cumulative impact on the final 
distance to HD 23642 does not exceed ~1 pc. The orbital 
solution is given in Table 2 and over plotted to the data 
in Figure 1. The parameters in common with the radial- 
velocities-only solution of Torres (2003) are very similar. 

The two stars are well within their Roche lobes and not 
distorted. The orbit appears relaxed to stationary condi- 
tions. In fact, the eccentricity is null within a small uncer- 
tainty and the rotational velocity of the two components, 
as measured on the high resolution spectra (37±1 for the 
primary and 31±2 km for the secondary), is that ex- 
pected for the rotation period to be synchronous with the 
orbital period. 

The barycentric radial velocity of HD 23642, 
y^=5.17±0.24 km s~^, is coincident with the mean ra- 
dial velocity of the Pleiades within the velocity disper- 
sion of its members. They are i?Vi;iustcr=5.19 km s^^ and 
(Trv=0.70 km s~^ according to Smith & Struve (1944), or 
i?V;iustcr=5.74±0.07 km s^^ and crRv=0.69±0.05 km 
following Narayanan & Gould (1999). 

The distance to HD 23642 resulting from our orbital 
solution is 132±2 pc (cf. Table 2). This result does not 



depend on the measured reddening for HD 23642. Rising 
it to i?B_\/= 0.035 average reddening for cluster members 
(e.g. Crawford & Perry 1976), would require an increase 
to 9910±145 K for the primary to fit observed colors (the 
larger error indicates a much poorer fit to literature 
photometry compared to i?B_y=0.012). The larger intrin- 
sic brightness compensates for the extra AAy=0.07 mag 
extinction and the resulting distance is 130.6±3.7pc. Even 
if the distance to HD 23642 derived in this paper con- 
cerns just one cluster member, nevertheless it is in excel- 
lent agreement with the results from the main-sequence 
fitting, moving cluster method and astrometric orbit of 
Atlas. These results agree so well that the issue of the 
distance to the Pleiades seems now addressed, and finally 
over. 
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